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Nuclei with large neutron-to-proton ratios have neutron skins, which manifest themselves in an 
excess of neutrons at distances greater than the radius of the proton distribution. In addition, 
some drip-line nuclei develop very extended halo structures. The neutron halo is a threshold effect; 
it appears when the valence neutrons occupy weakly bound orbits. In this study, nuclear skins 
and halos are analyzed within the self-consistent Skyrme-Hartree-Fock-Bogoliubov and relativistic 
Hartree-Bogoliubov theories for spherical shapes. It is demonstrated that skins, halos, and surface 
thickness can be analyzed in a model-independent way in terms of nucleonic density form factors. 
Such an analysis allows for defining a quantitative measure of the halo size. The systematic behavior 
of skins, halos, and surface thickness in even-even nuclei is discussed. 
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I. INTRODUCTION 

One of the main frontiers of nuclear science today is 
the physics of radioactive nuclear beams (RNB). Exper- 
iments with beams of unstable nuclei make it possible 
to look closely into many unexplored regions of the pe- 
riodic chart and many unexplored aspects of the nuclear 
many-body problem . 

Prospects for new physics, especially on the neutron- 
rich side of the beta-stability valley, have generated con- 
siderable excitement in the low-energy nuclear physics 
community. Neutron-rich nuclei offer an opportunity to 
study the wealth of phenomena associated with the close- 
ness of the particle threshold: particle emission (ion- 
ization to the continuum) and characteristic behavior 
of cross sections 1^,0], existence of soft collective modes 
and low-lying transition strength and dramatic 

changes in shell structure and various nuclear properties 
in the sub-threshold regime [p^-p^ . 

A very interesting aspect of nuclei far from stability 
is an increase in their radial dimension with decreasing 
particle separation energy [p7|-pl|. Extreme cases are 
halo nuclei - loosely bound few-body systems with about 
thrice more neutrons than protons. The halo region is a 
zone of weak binding in which quantum effects play a 
critical role in distributing nuclear density in regions not 
classically allowed. 

Halo nuclei, with their intricate topologies, are sym- 
bols of RNB physics. The very weak binding of the out- 
ermost neutrons leading to a rather good decoupling of 
halo from the core simplifies many aspects of underlying 
nuclear structure and reaction mechanism. Theoretically, 
the weak binding and corresponding closeness of the par- 
ticle continuum, together with the need for the explicit 
treatment of few-body dynamics, makes the subject of 



halos both extremely interesting and difhcult [Q. 

In the heavy neutron-rich nuclei, where the concept 
of mean field is better applicable, an a priori separation 
into core and halo nucleons seems less justified. However, 
the fact that there are far more neutrons than protons in 
these nuclei implies the existence of the neutron skin (i.e., 
an excess of neutrons at large distances). In addition, 
in neutron-rich weakly bound nuclei, one expects to see 
both the skin and the halo. 

There is no consensus in the literature on how to define 
and parametrize skins and halos. A quantity which is 
often employed to characterize the spatial extension of 
neutron density is the difference between neutron and 
proton root mean square (rms) radii: 



Ai?„p = {rl) - {rl). 



(1.1) 



In normal nuclei, this quantity is known to vary be- 
tween 0.1-0.2fm [p^-p5|, but it increases significantly in 
neutron-rich systems due to both skin and halo effects. 

Stimulated by recent experimental developments 
mainly in light nuclei, where some information on Ai?„p 
has been obtained [p6|-^ (see also the recent studies 
based on the giant dipole |^ and spin-dipole res- 
onance data and on antiprotonic levels |^,Q), many 
theoretical papers with a focus on radii of neutron and 




(For ear- 
) In most cases. 



proton density distributions appeared_ 
lier works, see papers quoted in Ref. 
theoretical studies were concerned with rms radii, an d 
the skin was usually discussed in terms of quantity ( |l.l| ) . 

Unfortunately, the second moment of nucleonic density 
(the rms radius) provides a very limited characterization 
of the nucleonic distribution. In particular, since the 
parameter Ai?„j, can be strongly influenced by weakly 
bound valence nucleons, i.e., by the shell structure, it is 
not able to properly describe the bulk radial behavior of 
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drip-line nuclei. A powerful tool that allows a more de- 
tailed description is the Helm model, introduced in the 
context of electron scattering experiments [|6-5l_ 
this model, the diffraction radius and surface thickness 
extracted from the density form factor are mainly sensi- 
tive to the nucleonic distribution in the surface region, 
and they are practically independent of shell fluctuations 
in the nuclear interior ||59|--|63|| . The robustness of the 
Helm model parameters, and their simple geometric in- 
terpretation, make this model a very attractive tool when 
characterizing density distributions. 

The main goal of this study is to apply the Helm model 
to nucleonic densities calculated in the self-consistent 
mean-field theory. In the first part of this paper, it is 
demonstrated that by analyzing the nucleonic form factor 
one is able to define, in a model-independent way, con- 
tributions to proton and neutron radii coming from skins 
and halos. In the second part, we perform systematic cal- 
culations of skins and halos in spherical even-even nuclei 
and discuss their dependence on the model employed. 

The material contained in this study is organized as 
follows. The analysis of nucleonic density based on the 
Helm model is outlined in Sec. ||. Section |l| discusses 
the details of Hartree-Fock-Bogoliubov (HFB) and rel- 
ativistic Hartree-Bogoliubov (RHB) models employed. 
The results of self-consistent calculations for diffraction 
radii, surface thickness, skins, and halos in spherical 
even-even nuclei are discussed in Sec. IV, together with 
the simple analysis based on the square-well model. Fi- 
nally, Sec. contains the main conclusions of this work. 



II. SPHERICAL HELM MODEL 

A key feature of the nucleonic density is the rms radius 



J (Pr p{r) 



(2.1) 



Further characteristics are best deduced from the corre- 
sponding form factor 



e"^'' p{r)d?r. 



(2.2) 



For the spherical density distribution /o(r), the form fac- 
tor F(g) is spherical and can be expressed in the standard 
way: 



F{q) = / jo{qr)p{r)r^dr. 



(2.3) 



There are various ways to characterize the basic pat- 
tern of the nucleonic density. A choice that is straight- 
forward and easy to use is provided by the Helm model 
p6|-p8[. Here, nucleonic density is approximated by a 
convolution of a sharp-surface density with radius Rq 
with the Gaussian profile, i.e. 



p(H)(r) = 



J d^r' fG{r-r')poe{Ro 



where 



fGir) 



(27r)3/2 



Ir'l), (2.4) 



(2.5) 



The radius in Rq in Eq. ( |2.4| ) is the diffraction (bo x equiv- 
alent) radius, and the folding width a in Eq. ( ^.5| ) models 
the surface thickness. The density pQ is given by 

hence the Helm density p*-^-* is normalized to the particle 
number N. The advantage of the Helm model is that 
folding becomes a simple product in Fourier space, thus 
yielding 



Roq 



(2.7) 



It is obvious that the first zero of F'-^'>{q) is uniquely 
related to the radius parameter Rq . The fit of this model 
parameter is thus trivial. We simply relate Rq to the first 
zero of the realistic form factor F{q), i.e., 

i?o = 4.49341/gi, (2.8) 

where qi is the first zero of F{q). This means that Rq 
can be deduced from the diffraction minimum and this is 
why it is called a diffraction radius (or the box-equivalent 
radius). The surface thickness parameter, ct, can be com- 
puted by comparing the values of microscopic and Helm 
form factors, F{qm) and F^^'^{q„i), at the first maximum 
qm of F{q), which gives 



SRhiiqM 



a^^^ln -"^^^;--^ . (2.9) 

qfn Roqm,F[qm) 

We have now at our disposal three key parameters that 
characterize the microscopic nucl eonic density: the rms 
radius -Rrms as defined in Eq. (2.1), the diffraction radius 
_Ro from Eq. (2.8), and the surface thickness a given by 
Eq. (2.9). The Helm model has only two independent 



parameters, and thus its rms radius can be expressed in 
terms of i?o and a: 



R^i^\l-{Rl + 5<j^ 



(2.10) 



Furthermore, it is more natural to discuss radii which 
pertain to a geometrical size of the nucleus and, there- 
fore, the prefactor in Eq. ( 2.10| ) is rather inconve- 
nient. We prefer to work with radii which we call here 
the geometric radii, defined as 



Rgcom — ^ g -Rrms J (2.11a) 
i?Hclm = = ^/(R^+5^). (2.11b) 

With this definition, the geometric radius becomes the 
box-equivalent radius in the limit of a small surface thick- 
ness. 
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III. MEAN-FIELD MODELS 

This section contains a very brief description of self- 
consistent models applied in this work. Since these mod- 
els are standard, our discussion is limited to basic defini- 
tions and references. 



A. Hartree-Fock-Bogoliubov model 

The HFB approach is a variational method which uses 
nonrelativistic independent-quasiparticle states as trial 
wave functions []6^ . An independent-quasiparticle state 
is defined as a vacuum of quasiparticle operators which 
are linear combinations of particle creation and annihila- 
tion operators. In this work, instead of using the matrix 
representation corresponding to a set of single-particle 
creation operators numbered by the discrete index, we 
use the spatial coordinate representation . This is 

particularly useful when discussing spatial properties of 
the variational wave functions and the coupling to the 
particle continuum |]6^ , p^ . 

In our HFB calculations, we employ the zero-range 
Skyrme interaction in the particle-hole channel. The 
total binding energy of a nucleus is obtained self- 
consistently from the energy functional |67|]: 



£ — £vi 



kin 



Sk 



£: 



Sk.ls 



pair 



£ 



CM, 



(3.1) 



where 5kin is the kinetic energy functional, £sk is the 
Skyrme functional, £sk,is is the spin-orbit functional, £c 
is the Coulomb energy (including the exchange term), 
£pa.ir is the pairing energy, and £cm is the center-of-mass 
correction. 

In this work, two Skyrme parametrizations are used: 
SkP and SLy4 |68|. Both of these selected forces per- 
form well concerning the total energies and radii. In par- 
ticular, both SkP and SLy4 parametrizations have been 
shown to reproduce long isotopic sequences p8| , |69[ . 

In the particle-particle channel, we use the SkP 
parametrization in the HFB/SkP variant and the 
density-independent volume delta interaction in the 
HFB/SLy4 variant. The strength of the delta force 
was adjusted according to the prescription given in 
Refs. II^,^. For the details of the calculations, we refer 
the reader to Refs. |6^,|Tf 



scalar cr-meson, the isoscalar vector w-meson, and the 
isovector vector p-meson. The model is based on the one- 
boson exchange description of the nucleon-nucleon inter- 
action. The starting point is the effective Lagrangian 
density |7l|,|| 

C — ^ («7 • d — m) ij) 

+ \{daf - U{a) - + \ral,j^ 

-^Rm.R^"^ + \mlp^ ~ Jf^.F^"^ - g^^ai^ 
-g^^p-f ■ ujip - gp^p-f ■ pTip - 6-07 • A^i^^V, (3-2) 



where 



U{a) 



r.920- 



(3.3) 



Vectors in isospin space are denoted by arrows. (Vectors 
in three-dimensional coordinate space are always indi- 
cated by bold-faced symbols.) The Dirac spinor ip rep- 
resents the nucleon with mass m, and mo-, m^, and rrip 
are the masses of the cr-meson, the w-meson, and the 
p-meson, respectively. The meson-nucleon coupling con- 
stants, go-, gun and gp, and unknown meson masses are 
parameters adjusted to fit nuclear matter data and some 
static properties of finite nuclei. U{a) denotes the non- 
linear cr self-interaction ||7| and , Rf"" , and F'"' are 
field tensors. 

For the purpose of the present study, we choose two 
RMF parameterizations: NL3 @ and NL-SH |Q. The 
force NL3 stems from a fit including exotic nuclei, neu- 
tron radii, and information on giant resonances. The NL- 
SH parametrization was fitted with a bias toward isotopic 
trends and it also uses information on neutron radii. 

The relativistic extension of the HFB theory was intro- 
duced in Ref. | |76[ | . In the Hartree approximation for the 
self-consistent mean field, one obtains the RHB equations 
which are solved sclf-consistently in coordinate space by 
discretization on the finite element mesh | |7^ ]. The spa- 
tial components, a;, ps, and A vanish due to the time- 
reversal symmetry. Because of charge conservation, only 
the third component of the isovector rho meson con- 
tributes. In the present investigation, the pairing in- 
teraction has been approximated by a phenomenological 
finite-range Gogny force with the DIS parameter set [ [78[ . 
This force has been adjusted to the pairing properties of 
finite nuclei all over the periodic table. 



B. Relativistic mean- field model 



IV. RESULTS AND DISCUSSION 



Relativistic mean-field theory has been proved to be 
a powerful tool in describing various aspects of nuclear 
structure Q . The model explicitly includes mesonic de- 
grees of freedom and describes the nucleons as Dirac par- 
ticles. Nucleons interact in a relativistic covariant man- 
ner through the exchange of virtual mesons: the isoscalar 



A. Skins and halos in spherical heavy nuclei 

The Helm-model characteristics of calculated density 
distributions are obtained from the microscopic form fac- 
tors (2.3). Figure |l| shows the neutron densities for 
I20,i5cfrr7iisn calculated in the RHB/NL3 model and the 
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corresponding form factors. The positions of the first 
and second zero of the form factor (indicated by ar- 
rows) decrease gradually with neutron number refl ectin g 
the steady increase of the neutron radius [see Eq. (2.8)]. 
The zeros of F{q) are regularly spaced and the ratio of 
(72/91 is very close to the ratio of the flrst two zeros of 
the spherical Bessel function ji. It is also seen that, 
in the considered range of q, the envelope of q\F{q)\ is 
practically constant Q . All of these observations con- 
firm that in the region of \ow-q values shown in Fig. ||, 
the "model-independent" analysis of theoretical density 
distributions, according to Ref. |60|, can safely be per- 
formed. 

Our analysis of neutron and proton radii in the Sn 
isotopes is summarized in Fig. y. The most interesting 
observation is that for the isot opes w ith A^<82, the neu- 
tron geometric radiu s i^eco m ( ^.lla ) is very close to the 
Helm radius i?Hcim (2.11b). On the other hand, for nu- 
clei heavier than -'^■^^Sn, the former is appreciably greater 
than the latter. This behavior suggests that the differ- 
ence between i?goom and -Rncim is related to the size of the 
neutron separation energy. Indeed, for iV<132, the two- 
neutron separation energy is S'2n>12MeV, and it drops 
to a few MeV around iV=100. Due to the weaker binding, 
the neutron distributions in the very heavy tin isotopes 
have larger spatial extensions, and this increases i?gcom 



dramatically due to the weig ht in Eq. (|2.l|). On the 



other hand, the form factor at intermediate values of q is 
almost independent of the asymptotic tail of the density 
distribution. Therefore, the radius parameters deduced 
from the form factor, Rq and i?Hcim, show a less dramatic 
growth. 

Guided by this observation, we introduce the halo pa- 
rameter as the difference 



SRhalo — R. 



gcom 



-R 



Holm- 



Such a halo parameter is indicated in Fig. ^, where it 
shows the size of the neutron halo in neutron-rich tin 
isotopes. It should be noted that a halo may also be 
defined through the higher radial moments, e.g., (r*). 
We have checked, however, that other definitions d o no t 
have any advantage over the simplest prescription ( [4.l|) . 

In contrast to the neutron halo parameter i5i?haio(n), 
for protons the value of (5i?haio(p) turns out to stay very 
close to zero, i.e.. 



R^ 



1 (P) ~ -RHclin(p), 



(4.2) 



such that one cannot easily resolve the difference of 
-Rgcom(p) and i?Hcim(p) in the plot. This reflects the fact 
that protons are always very well localized in the nu- 
clear interior by the Coulomb barrier, and they are very 
well bound. (The two-proton separation energy increases 
from S'2pw4MeV in i^^Sn to «28MeV in "OSn.) 

Figure || shows the calculated neutron densities for the 
tin isotopes multiplied by r**. (The area under pr'' is 
proportional to (r^).) It is seen that the large value of 
-Rgeom in HFB/SLy4 can be attributed to the presence of 



a hump in pr"* at rwlS fm. In the particular represen- 
tation of Fig. ^, it is possible, at least in the HFB/SLy4 
model, to see rather clearly the decoupling of nuclear den- 
sity into the "core" and "halo" parts. This decoupling 
is much weaker in the HFB/SkP model, and it is almost 
invisible in RHB/NL3. 

While the halo is a property of neutrons or protons, the 
neutron skin depends on the difference between neutron 
and proton radii, and thus it is more difficult to quantify. 
Indeed, since several deflnitions of a radius have been 
employed in this work, one can introduce various param- 
eters reflecting the neutron- proton radius difference, e.g.. 



Ai?o 

Ai?Helm 



AJ?geom — Rgeom (n) 



i?o(n)-i?o(p), 

-RHelm(n) — -RHclin(p), 
Rgcom (p) • 



(4.3a) 
(4.3b) 
(4.3c) 



These three definitions are displayed in Fig. ^ for the Sn 
isotopes calculated in the IIFB/SLy4 model. 

According to the discussion above, the difference of ge- 
ometric radii, Ai?gcom, contains a contribution from halo 
effects; hence it is not appropriate to define the skin. 
The differences ARq and Ai?Hcim both smoothly increase 
with neutron number, with Ai?Hcim being always greater 
than Ai?o, due to the contribution from the surface thick- 
ness. In principle, both definitions could be used to char- 
acterize the skin. However, due to the smallness of the 
proton halo ([l.2|), one simply has 



for 



Ai?goom ~ '^^halo(n) + SRskin 



in ; 



(4.4) 



(4.5) 



(4-1) i.e., Eq. (4.4) gives an additive decomposition into the 



contributions to Ai?geom coming from the weak binding 
(halo part) and representing the size effect (skin part) 
(see Fig. Therefore, for the present purpose, we prefer 
Ai?Hcim as a measure of the skin. Figure | nicely shows 
that the neutron halo effect in the Sn isotopes is predicted 
to show up just above N=82, and it increases gradually 
with TV reaching in the HFB/SLy4 calculations the value 
of (5i?haio(n)«0.65Ai?Hoim near the two-neutron drip line. 

The results of calculations for the Ni isotopes are 
shown in Fig. |^. Here, the neutron skin quickly in- 
creases above the doubly magic nucleus ^^Ni, i.e., above 
the N—50 gap. A simpler pattern is seen for the Pb iso- 
topes (see Fig. ||): the neutron halo develops for N>126. 
In all cases, <5i?haio calculated with SLy4 is systemat- 
ically greater than that in HFB/SkP, RHB/NL3, and 
RHB/NLSH. 

In order to understand these results, we show in Figs. |^ 
and ^ the two-neutron separation energies, S2n, for the 
neutron-rich Ni, Sn, and Pb isotopes. Systematically, the 
HFB/SLy4 model predicts the lowest separation energy. 
When approaching the neutron drip line, both RHB ap- 
proaches yield considerably larger neutron binding than 
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Skyrme-HFB calculations. This result is consistent with 
the model dependence of i5i?haio- Indeed the neutron halo 
parameter seems to be correlated with the neutron sepa- 
ration energy. Tha t is, (5i?haio increases with decreasing 
82,1 (see Sec. IVC for more discussion concerning this 
point). 



The surface thickness cr (2.9) shows characteristic de- 
pendence on particle number (see Fig. ||). Namely, cr„ 
increases with N on the average, but it shows local min- 
ima around magic numbers. This local decrease in a can 
be attributed to its sensitivity to pairing correlations . 
Indeed, static pairing correlations in magic nuclei vanish, 
the Fermi surface becomes less diffused, and the surface 
thickness is reduced. When approaching the spherical 
neutron drip line, cr„ behaves fairly smoothly; because 
it is determined from the formfactor at larger q (i.e., it 
seems to be rather insensitive to the asymptotic behav- 
ior of nucleonic density at large distances). The proton 
surface thickness behaves fairly constant as a function of 
TV, although it also exhibits the local decrease at magic 
neutron numbers as a result of self-consistency. 

Except for the very neutron-rich nuclei, the RHB mod- 
els yield tr-values which are lower than in the Skyrme- 
HFB calculations. This effect is particularly clear for 
CTp, which is not affected by the variations in the pairing 
field. In addition, in all cases crp(NLSH)<(Tp(NL3) and 
CTp(SLy4)<crp(SkP). (For further discussion, we refer the 
reader to Ref. @.) 

The difference (T„-crp exhibits very weak shell effects. 
It gradually increases from about 0.2 fm around the beta 
stability line to about 0.5 fm near the neutron drip line. 
Interestingly, as discussed in Ref. the difference be- 
tween neutron and proton radii also depends very weakly 
on shell effects. 



B. Square- well potential analysis 

This section contains some general arguments regard- 
ing the concept of diffraction radius in two-body halo 
systems. Our discussion is based on the spherical finite 
square-well (SW) potential used in Ref. |ll7|l to illustrate 
some generic aspects of halos. (See Ref. |p2|| for the ex- 
tension to the deformed case.) 

The advantage of this simple model is that by changing 
the well depth, one can vary the position of bound single- 
particle halo orbitals and, therefore, study the properties 
of diffraction radii and surface thickness very close to the 
e=0 threshold. In our calculations we assume that the 
square- well potential radius is i?sw=7fm and the system 
consists of 70 particles. The potential depth is varied to 
tune the energy of the last bound nuclcon. The halo 
structure is represented by two neutrons in the 3si/2 or- 
bital, while the core can be associated with the remaining 
68 particles occupying well-bound states. 

With the binding energy of the 3si/2 orbital approach- 
ing zero, the halo develops. This is illustrated in Fig. nO, 



which shows the total and core densities for three values 
of the binding energy of the 3si/2 halo orbital: -5MeV, 
-100 keV, and -IkeV. The presence of the halo is clearly 
seen at larger distances, r>8fm. 

The form factors of the total and core densities at 
e(3si/2)=--5MeV and -IkeV are displayed in Fig. |ll]. As 
a consequence of the uncertainty principle, in the case of 
a very weak binding the form factor of the halo wave func- 
tion (shown in the insert) corresponds to a very narrow 
momentum distribution, and it contributes very little to 
the total form factor. Consequently, the first zero of the 
form factor, hence the diffraction radius, is very weakly 
influenced by the presence of the halo. This is not true 
in the case of a large binding where the valence orbital 
docs not have the halo character, and its form factor is 
significantly greater from zero in the region of qi . 

Figure ^ displays the radii calculated as a function 
of the 3s 1/2 binding energy. Due to the halo character 
of the valence orbital, with the energy of the 3s 1/2 state 
approaching zero, the total geometric radius diverges as 
(— e)~^ At the same time, the geometric radius of 

the core, as well as the Helm radii for the total system 
(-RHoim.totai) and the core (i?Hoim,core), very weakly de- 
pend on e. The effect of decreased binding on the core is 
measured by the difference i?geom, core— -Rncim, core, which 
gives the core contribution to (5i?haio- As expected, in 
the limit of weak binding, the halo parameter is almost 
entirely determined by the asymptotic behavior of the 
3s 1/2 wave function. It is also seen that the difference 
between i?Heim,totai and -RHeim.core is very smaU. 



C. Pairing anti-halo effect 

In contrast to light nuclei where the halo can be as- 
sociated with very few weakly bound neutrons that are 
practically decoupled from the rest of the system, it is 
difficult to separate halo structures in heavier systems 
in which all the nucleons (including the valence weakly 
bound neutrons) move in one self-consistent field. An 
additional difference and complication is caused by the 
presence of strong pairing correlations in heavy open- 
shell nuclei. As found in Ref. and discussed below, 
pairing strongly modifies the extreme si ngle-p article pic- 
ture of halo structures presented in Sec. IV B . 

Consider, e.g., the valence neutron moving in a mean- 
field potential. Due to the fact that the nuclear mean 
field vanishes at large distances, the standard asymptotic 
behavior of the neutron density is, in the absence of pair- 
ing correlations, given by 



p{r) 



exp (— 2Kr) 



where 



l2m{-e) 



(4.6) 



(4.7) 



5 



with e being the single-particle energy of the least bound 
neutron. In the presence of pairing, the constant k is 
831 , and for the even neutron numbers 



different 
reads 



5N66 16 



(4.8) 



where Emin is the lowest quasi-particle energy and A is 
the Fermi energy. 

In the extreme single-particle picture (no pairing), the 
halo structure may develop when e — > (/c — > 0) [0. 
However, according to expression (4.8), in the limit of 
vanishing binding (A —^ 0) the constant n does not vanish 
and reads 



2mA 



(4.9) 



where A is the pairing gap of the lowest quasiparticle. 
Consequently, in the presence of pairing correlations, k 
is never small, and a huge halo, as it is seen in light nuclei, 
cannot develop {pairing anti-halo effect oi Kef. ||8^). 

In order to confirm the influence of pairing and weak 
binding on SRhaio, we performed spherical IIFB/SLy4 
calculations near the two-neutron and two-proton drip 
lines. Figure [l^ shows the neutron halo parameters, neu- 
tron Fermi energies, and neutron pairing gaps calculated 
in the IIFB/SLy4 model for the two- neutron drip-line 
even-even nuclei (which are the heaviest even-even iso- 
topes that are still predicted to be two- neutron bound). 
First, we note that A„ does not vanish near the two- 
neutron drip line. This phenomenon has been found and 
discussed in detail in Refs. [IbF 



831, and it was at- 
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tributed to the strong coupling to the neutron contin- 
uum in the pairing channel. Second, the pairing gap 
shows some shell fluctuations: the minima in A„ ap- 
pear at N=86, 130, and 192, i.e., just above the neutron 
magic gaps. On the average, however, A„ stays between 
'--^l.S MeV in light nuclei and ^^1.2 MeV in the heaviest el- 



ements. As a result, the exponent (4.9) is always sizable, 
and (5i?haio (ii) does not exceed 1 fm in heavy even-even 
nuclei. 

The pattern of (5i?haio(n) seen in Fig. |l^ is nicely corre- 
lated with the behavior of A„ . Namely, the neutron halo 
parameter increases when the Fermi energy approaches 
zero. It is to be noted, however, that there is no clear 
correlation between the magnitude of i5i?haio(ii) and the 
appearance of \ow-£ (s and p) states at the Fermi energy 
I prf , as one would expect from an extreme single-particle 
picture of Sec. IV B . It seems that the pairing anti-halo 
effect is far more important than the influence of the 
centrifugal barrier, cf. discussion in Ref. We made 

an attempt to find a phenomenological expression that 
would express i5i?haio(n) in terms of (CA„ — A„)'' (^,7y 
being free parameters). Unfortunately, we were not able 
to obtain a unique fit for all neutron-weak nuclei at once, 
although some correlation between these two quantities 
exists. 



The insert in Fig. |l^ shows the proton halo parameter 
in the least bound even-even isotones near the two-proton 
drip line. As expected, due to the confining effect of the 
Coulomb potential, the proton halo is very small - of the 
order of 0.02 fm. It is only in the very light sd nuclei that 
(5i?haio(p) can exceed 0.1 fm. Interestingly, there is also 
some increase in the proton halo in the superheavy nuclei 
with Z^120, A^~172, which, in some spherical calcula- 
tions, show bubble-like structures In this context, 
it should be emphasized again that calculations shown 
in Fig. |l^ are spherical, and some modifications due to 
deformation are expected; in particular, the superheavy 
nucleus with Z=120 and iV=172 is not expected to be 
spherical in the HFB/SLy4 model 



D. Global behavior of halos and skins in spherical 
even-even nuclei 

In order to study the systematic behavior of the 
spherical-shape density distributions, we performed sys- 
tematic calculations in the spherical IIFB/SLy4 and 
HFB/SkP models for all even-even nuclei predicted to 
be stable with respect to the two-nucleon emission, 
i.e., for all even-even nuclei with positive two-neutron 
and two-proton separation energies, S2n=B(N, Z) — 
B{N-2, Z)>0 and S2p^B{N, Z) - B{N, Z-2)>0. 

The results for the neutron halos are shown in Figs. |lj 
and|l|for the HFB/SLy4 and HFB/SkP models, respec- 
tively. Several features seen in these systeniatics are note- 
worthy. Firstly, for most nuclei 5i?haio(n) is very small. 
Only in the immediate vicinity of the two-neutron drip 
line is a rapid increase in the halo parameter seen. As 
discussed above, while the halo effect is rather strong for 
the SLy4 force, the HFB/SkP model predicts very few 
candidates for a halo. 

A second interesting aspect is a weak dependence of 
the halo parameter on shell effects. Contrary to the rms 
radii which show a significant reduction around spher- 
ical magic gaps p5[ |, the variations of (5_Rhaio(ii) around 
magic gaps are much weaker. This is easy to understand. 
In well-bound nuclei where the shell effects are very pro- 
nounced, the halo parameter is dramatically reduced due 
to weak binding. On the other hand, in neutron drip-line 
nuclei, shell effects are significantly weakened (reduced 
magic gaps, strong pairing correlations); hence their in- 
fluence on radii is less significant. The pattern shown 
in Figs. |lj and |l^ basically reflects the behavior of the 
two-neutron separation energy around the two-neutron 
drip line, and is qualitatively similar to that for Ai?np 
discussed in Rcf. [p5|. 

As shown in Figr]l3|, the proton halo parameter is much 
smaller than that for the neutron halo. The inset in 
Fig. |lj shows (5i?haio(p) for very light nuclei. The largest 
proton halos, (5_Rhaio(p)'~-"0.15 fm, can be fo und a round 
20Mg and ^''Si. For more discussion, see Sec. fVE. 
shows the 



Figure shows the neutron skins calculated in the 
HFB/SLy4 model. One sees that the skin grows steadily 



6 



in a direction orthogonal to the valley of stability. The 
weak mass dependence and a nearly linear trend with 
the neutron excess N—Z suggests that <5i?skin reflects 
the bulk size properties of neutrons and protons. As dis- 
cussed in Ref. , the isovcctor dependence of the neu- 
tron skin is governed by a balance between the volume 
pl| , p2t and surface symmetry energy coefficients. Within 
the present sample, this is confirmed by the fact that the 
HFB/SkP results for the neutron skin are indeed very 
similar, and it is noted that SkP and SLy4 do have a very 
similar symmetry energy coefficient (32 MeV in SLy4 and 
30MeV in SkP). 

Last but not least, it is worth inspecting the global 
trends of the surface thickness CTp n. Figure |l^ shows the 
neutron surface thickness. It s pattern shares one feature 
with the neutron halo (Fig. p^ , namely dn is particularly 
large near the neutron drip line. However, dn displays a 
much richer structure all over the periodic table with 
maxima far from closed shells. The proton surface thick- 
ness is shown in Fig. Compared to (T„, the global 
behavior of CTp is different. The variations with N—Z are 
less systematic and, in some cases, ap decreases when ap- 
proaching the two-proton drip line. As in neutrons, the 
proton surface thickness is reduced around magic gaps. 
As discussed earlier, a-p is generally much smaller than 

(Tn- 

It is interesting to note the presence of an island of par- 
ticularly small neutron and proton skins near the proton 
drip line in the region of superheavy nuclei with A^~172. 
This is probably related to the pronounced dip of the 
spherical distribution near the nuclear center which ap- 
pears for these nuclei |^ , ^ . For the protons, there ex- 
ists a further island of small surface thickness for super- 
heav y elem ents with iV>200. However, as discussed in 
Sec. [VC, the presence of bubble- like structures in this 



region may be an artifact of the assumption of spherical 
symmetry jo^ . 

Finally, the neutron-proton difference of the surface 
thickness is shown in Fig. ^ It displays a mix of steady 
growth with A as well as with N — Z which is just the 
sum of the different trends seen for protons and neutrons 
separately. Like the diffraction radii and rms (or geo- 
metric) radii, the shell effects which are present in both 
observables are somewhat suppressed in the differential 
quantity. 



model, such a situation was found for the neutron-rich 
Ne isotopes. According to this work, the neutron Fermi 
energy in the Ne isotopes with A^>20 stays very close 
to zero, stabilized by the presence of three close-lying 
single-particle canonical orbitals, 2^3/2, 2pi/2, and 1/7/2- 

The low-^ shells, 2^3/2 and 2pi/2, are good candidates 
for halo orbitals. It is worth noting that their canoni- 
cal HFB energies stay very close in energy. This sug- 
gests that the wave functions of the 2^3/2 and 2pi/2 are 
weakly influenced by the spin-orbit interaction (a situ- 
ation that is characteristic of halo states) ||9^. This is 
nicely demonstrated in Fig. ^ which shows the form fac- 
tors of single-neutron canonical states lp3/2j 1c^3/2j 2p3/2, 
and 2pi/2 in the drip- line nucleus ^^Ne. The form fac- 
tors of the 2p3/2 and 2pi/2 orbitals are very similar which 
confirms the negligible effect of the spin-orbit interaction 
on these states. The narrow momentum distribution of 
the 2p orbitals is indicative of weak binding. In contrast, 
wider form factors of well-bound IJ33/2 and 1^3/2 orbitals 
reflect the fact that their wave functions are better local- 
ized inside the nuclear volume. 

The neutron distribution form factors in ^"'■^'^Ne ob- 
tained in RHB/NL3 are shown in Fig. |l|. In ^ONe, the 
contribution from the valence 1^5/2 neutrons has been 
singled out, and it is seen that the influence of the va- 
lence orbits on the diffraction radius is strong. On the 
other hand, the effect of the 2p valence orbitals on Rq 
in •^"'Ne is small. The reason for this is twofold. Firstly, 



in accordance with the discussion from Sec. IV B, the 2p 
form factor is narrow and it mainly contributes around 
q—Q. Secondly, there are less than two neutrons in the 2p 
shell. Hence the behavior of the total form factor in '^''Ne 
(the ratio of the number of neutrons in valence orbits to 
that in the core is ~0.08) is primarily governed by the 
core neutrons. (In ^"Ne the valence/core ratio is 0.25.) 

The above situation discussed for '^^Ne is, in fact, typ- 
ical for all weakly bound neutron-rich nuclei. The weak 
binding of valence orbits and the fact that they are oc- 
cupied by very few particles makes the diffraction radius 
weakly dependent on halo structures. It is interesting 
to see that the simple argument presented in Sec. IV B 



works in a microscopic case, in spite of the fact that pair- 
ing modifies the naive single-particle picture to some ex- 
tent. 



E. RHB calculations of halos in the Ne isotopes 



In contrast to the simple model of Sec. |lVB or analysis 
of Ref. | |8^ , in microscopic calculations the binding en- 
ergy of a single particle (or quasi-particle) orbital is not a 
free parameter but is obtained self-consistently from the 
realistic Hamiltonian. Hence it is difficult to find a case 
where a low-£ orbital (i.e., a potential candidate for halo) 
appears very close to the threshold. Here we discuss the 
case studied in Ref. |M, where, based on the RHB/NL3 



F. Charge halo parameter in stable nuclei 

Elastic electron scattering has provided a world of well- 
evaluated data on nuclear charge distributions; see, e.g., 
Ref. [|60y96[ . It is nicely corroborated by the very pre- 
cisely measured root- mean-square charge radii pTf . This 
offers a possibility to deduce the experimental charge halo 
parameter. 



<5-Rhalo(ch) = i?geom(ch) - i?Helm(ch), 



(4.10) 



in selected cases. In order to have a most complete and 
up-to-date supply of data on charge radii and surface 
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thickness, we have recurred to the data base of Ref. [ p8[ , 
which is a compilation of results analyzed as explained 
in (§]. 

The charge form factor Fq (q) is composed of the form 
factors of the proton and neutron distribution multi- 
plied by the corresponding intrinsic nuclcon form factors. 
A similar contribution from the magnetic form factors 
is added. Finally, the centre-of-mass correction is per- 
formed. The nucleon form factors are taken from elec- 
tron scattering on the proton and the deuteron [p9,100 
and parametrize d in terms of the Sachs form factors as 
outlined in Ref. |101]. A detailed description of the pro- 
cedure used to det ermi ne the charge radii can be found 
in Appendix 2 of |102|. From the charge form factor, 
we deduce the rms radius, the diffraction radius, and the 
surface thickness in the standard manner. 

Figure shows the experimental values of i5i?haio(ch) 
for selected nuclei (doubly magic ^'^O, ^^•'^^Ca., ^*Ni, 



Pb, semi-magic ^-^Cr, ^^Pe, ««Sr, «"Zr, ^^Mo 



116,124 



Sn, 



204,206p^^ and some open-shell nuclei, including the well- 
deformed Cr and Sm isotopes). They are displayed to- 
gether with predictions of spherical Skyrme HFB calcu- 
lations. 

It is to be noted that the charge halo is a very sen- 
sitive obser vable because it stems from subtracting two 
large radii ( 4.10 ). The experimental error on (5-Rhaio(ch) 
is at least as large as the largest error on radii and surface 
thickness. This leads to a conservative uncertainty on the 
data, ~0.03fm. As expected from our results on proton 
halos, the charge halos are all very small. A notable ex- 
ception is ^^O where 5i?haio(ch) is 0.13 fm. Also, our cal- 
culations are expected to slightly underestimate charge 
halos in some open-shell nuclei (e.g., ^^^'-^^"'Sm) due to 
possible contributions from deformation effects. Consid- 
ering the above, it is very satisfying to see that our HFB 
results are generally close to the experimental points, in 
fact staying within the experimental uncertainty in most 
cases. Since proton halos are close to the charge ha- 
los in all calculations, this nice agreement, together with 
the systematic behavior of proton halos discussed above, 
make us conclude that the pronounced proton halos do 
not exist. 

At second glance, one is tempted to spot shell effects 
and isotopic trends when looking at the fluctuation of the 
charge halos in Fig. ^ However, these variations stay 
within the experimental uncertainties and cannot serve 
for a deeper analysis. The charge halos as such are too 
small and, moreover, the regime of stable nuclei docs not 
supply enough variation for that. 



V. CONCLUSIONS 

This work contains the theoretical analysis of neu- 
tron and proton skins, halos, and surface thickness ob- 
tained within the spherical self-consistent mean-field the- 
ory. The main goal was to describe spatial character- 



istics of nucleonic densities of nuclei far from stability, 
where the closeness of the particle continuum qualita- 
tively changes the physical situation. 

The Helm model analysis presented in this work al- 
lows for an unambiguous determination of halos and skins 
from nucleonic density distributions. It has been shown 
that the halo parameter, defined as the difference be- 
tween the geometric radius (a rescaled rms radius) and 
the Helm radius, is small in well-bound nuclei, and for 
neutrons it becomes enhanced for heavy exotic systems 
with low neutron separation energies. However, unlike 
in the neutron-rich few-body systems, our calculations 
do not predict giant neutron halos in medium-mass and 
heavy nuclei. This is because strong pairing correla- 
tions effectively reduce the impact of weak binding on the 
asymptotic behavior of the single-particle density (pair- 
ing anti-halo effect [s^]). 

No significant proton halo has been found when ap- 
proaching the proton drip line. A moderate effect (less 
than 0.2 fm) is predicted for some light nuclei, but it can 
be practically neglected for heavier systems. The exper- 
imental values of charge halos for stable nuclei, of the 
order of 0.02-0.04 fm, are perfectly consistent with the 
mean-field predictions. 

The neutron skin, defined as a difference of neutron 
and proton Helm radii, shows a smooth gradual depen- 
dence on the neutron excess and is extremely weakly af- 
fected by shell effects. This is consistent with the results 
of a previous study | p5[ | where a very weak shell depen- 
dence of Ai?„p was found. 

On average, the neutron surface thickness increases 
with neutron number, but it is locally reduced around 
magic numbers, thanks to reduced pairing. On the other 
hand, proton surface thickness depends to a lesser degree 
on proton number; it rather tends to follow the trend dic- 
tated by f7„. As a result, the difference cr„ — ap shows a 
reduced dependence on shell effects. A very interesting 
situation is predicted for the superheavy A^=172 isotones 
where the proton surface thickness is actually reduced 
with increasing proton number. 

Theoretically, the analysis based on density form fac- 
tors is very simple and physically elegant. Unfortunately, 
experimental information on nucleonic densities is cur- 
rently limited to charge densities in some stable nuclei, 
and almost nothing is known on neutron density distri- 
butions. Starodubsky and Hintz |103| made an attempt 
to deduce neutron densities in 206,2n77208p|-, _ ^ model- 
dependent way - from elastic proton scattering at inter- 
mediate energies, and this marks the state-of-the-art. An 
exciting new avenue is a prospect for direct measurements 
of the neutron density form factors from the asymmetry 
in the parity-violating elastic polarized electron scatter- 
ing |104-10£]. Table shows the Helm model analysis of 
neutron densities in ''"Ca and ^°^Pb calculated in our self- 
consistent models. Theoretical predictions for diffraction 
radii are rather robust, with the differences between val- 
ues of Rq obtained in different models being below 0.5% 
for ^''Ca and below 3.5% for ^°®Pb. Interestingly, due to 
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a compensation effect between Rq and cr, the geometric 
neutron radius in ^"^Pb is rather similar in aU models, 



i?2 



i7.32 fm. We hope that our results will stimulate 



future experimental studies of neutron distributions in 
nuclei. 
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TABLE I. Characteristics of neutron distributions in ''"Ca 
and ^'^*Pb: diffraction radius _Ro, surface thickness a, and 
geometric radius 7?geom (all in fm), obtained in the HFB and 
RMF models employed in this work. 



nucleus 




SLy4 


SkP 


NL3 


NLSH 




Ro 


3.827 


3.844 


3.844 


3.841 


«Ca 


a 


0.905 


0.923 


0.845 


0.793 




-Rgcom 


4.353 


4.388 


4.296 


4.274 




Ro 


6.870 


6.849 


7.076 


7.075 


208pb 


On 


1.022 


1.033 


0.971 


0.929 




^gcom 


7.252 


7.244 


7.409 


7.374 
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FIG. 2. Neutron -Rgcom(n) and _RHcim(n), and pro- 
ton /?geom(p) radii for the Sn isotopes calculated in the 
HFB/SLy4, HFB/SkP, RHB/NLSH, and RHB/NL3 mod- 
els. The proton Helm radius -RHoim(p) is also shown in the 
RHB/NL3 variant (dotted line); it is very close to -Rgoom(p). 



FIG. 4. The neutron-proton radius differences (4.3) for the 
even-even Sn isotopes calculated in the HFB/SLy4 model. 



12 



C<5 



Ni 



HFB/SLy4 




HFB/SkP 



30 40 50 80 
Neutron Number 





F1HB/NL3 

■30' ' ' '40 ' ' '50' ' ' 'eo' ' ' 
Neutron Number 



FIG. 5. Same as in Fig. ^ except for the Ni isotopes. 
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FIG. 6. Neutron and proton radii for the Pb isotopes cal- 
culated in the HFB/SLy4 and HFB/SkP models. 
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FIG. 7. Two-neutron separation energies for the neu- 
tron-rich Ni (top) and Sn (bottom) isotopes calculated in the 
HFB/SLy4, HFB/SkP, RHB/NLSH, and RHB/NL3 models. 
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FIG. 8. Two-neutron separation energies for the neu- 
tron-rich Pb isotopes calculated in the HFB/SLy4 and 
HFB/SkP models. 
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FIG. 9. Neutron (top) and proton (middle) surface thick- 
ness coefficients for the Ni (left) and Sn (right) isotopes 
calculated in the RHB/NL3, RHB/NLSH, HFB/SkP, and 
HFB/SLy4 models. The difference, a„—ap, is shown in the 
bottom panels. 
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FIG. 10. Total (dotted line) and core (solid line) densities 
for the system of 70 particles moving in the the finite spherical 
square- well potential with the radius i?sw=7 fm for the three 
values of the binding energy of the 3si/2 halo orbital: -5 MeV 
(top), -100 keV (middle), and -IkeV (bottom). 



FIG. 11. Form factors of the total (dotted line) and core 
(solid line) densities for the system of 70 particles moving 
in the finite spherical square-well potential with the radius 
i?sw=7fm for the two values of the binding energy of the 
3si/2 halo orbital: —5 MeV (top) and — IkeV (bottom). The 
form factor of the halo wave function is shown in the inset. 
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FIG. 12. Dependence of the geometric radii, namely the 



total, R, 



■gcom, total , 



and core, R) 



geom,core 



radii and of the corre- 



sponding Helm radii, -Rncim.totai and 7?Heim,core, on the bind- 
ing energy of the 3si/2 halo orbital of the finite square well 
with radius 7?sw=7fm. The total number of particles is 
N=70. The core consists of 68 particles occupying all the 
single-particle orbitals below 3si/2. 
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FIG. 13. Neutron halo parameters (top), neutron Fermi energies (middle), and neutron pairing gaps (bottom) ealculated in 
the HFB/SLy4 model for the two-neutron drip-line even-even nuclei (i.e., the heaviest even-even isotopes which are predicted 
to be two- neutron bound). The proton halo parameters for the two-proton drip-line even-even nuclei are shown in the inset. 
Note the correlation between 5i?haio and An, marked by the vertical dashed lines; whenever An approaches zero, SRhaio tends 
to increase. 
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FIG. 14. Neutron halo parameters ( [l.l[ ) calculated in the spherical HFB/SLy4 model for two-particle stable even-even nuclei. 
The inset shows proton halo parameters for very light nuclei. 
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FIG. 15. Neutron halo parameters (|4.l|) calculated in the spherical HFB/SkP model for two-particle stable even-even nuclei. 
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FIG. 16. Skin parameters (4.5) calculated in the spherical HFB/SkP model for two-particle stable even-even nuclei. 
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FIG. 17. Neutron surface thickness calculated in the spherical HFB/SkP model for two-particle stable even-even nuclei. 
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FIG. 18. Same as in Fig. except for the proton surface thickness. 
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FIG. 19. Same as in Fig. [LTI except for the difference (j„ — ap. 
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FIG. 20. Form factors of canonical RHB/NL3 sin- 
gle-neutron orbitals IP3/2, 1^3/2, 2p3/2, and 2pi/2 in ''^Ne. 
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FIG. 21. Total (solid line), core (dashed line), and valence 
(dash-dotted line) neutron form factors for ^"Ne (top) and 
^''Ne (bottom). The valence space is here defined to be given 
by two lds/2 neutrons in ^"Ne and all the occupied 2p states 
in ^^Ne. The form factors at q=0 are normalized to the cor- 
responding neutron numbers (e.g., _Ftot(0)=A^). 
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FIG. 22. Comparison of experimental and theoreti- 
cal charge halos for a selection of stable nuclei. The spher- 
ical theoretical results are produced with the two Skyrme 
parametrizations SLy4 and SkP used throughout this paper. 
The estimated experimental errors are about 0.03 fm. 
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